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Abstract. Isothermal rebreathing has been proposed as an experimental technique 
for estimating the alveolar levels of hydrophilic volatile organic compounds (VOCs) 
in exhaled breath. Using the prototypic test compounds acetone and methanol we 
demonstrate that the end-tidal breath profiles of such substances during isothermal 
rebreathing show a characteristic increase that contradicts the conventional pulmonary 
inert gas elimination theory due to Farhi. On the other hand, these profiles can reliably 
be captured by virtue of a previously developed mathematical model for the general 
exhalation kinetics of highly soluble, blood-borne VOCs, which explicitly takes into 
account airway gas exchange as major determinant of the observable breath output. 
This model allows for a mechanistic analysis of various rebreathing protocols suggested 
in the literature. In particular, it predicts that the end-exhaled levels of acetone 
and methanol measured during free tidal breathing will underestimate the underlying 
alveolar concentration by factor of up to 1.5. Moreover, it clarifies the discrepancies 
between in vitro and in vivo blood-breath ratios of hydrophilic VOCs and yields further 
quantitative insights into the physiological components of isothermal rebreathing and 
highly soluble gas exchange in general. 
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1. Introduction 

Over the past decade, much advance has been made in the attempt to extract the 
diagnostic and metabolic information encapsulated in a number of endogenous volatile 
organic compounds (VOCs) appearing in human exhaled breath [Ifj3]. However, the 
quantitative understanding of the relationships between breath concentrations of such 
trace gases and their underlying systemic levels clearly lags behind the enormous 
analytical progress in breath gas analysis. In particular, formal means for evaluating 
the predictive power of various breath sampling regimes are still lacking. 

Recently, several efforts have been undertaken to complement VOC measurements 
with adequate physical models mapping substance-specific distribution mechanisms in 
the pulmonary tract as well as in the body tissues. Some major breath VOCs have 
already been investigated in this form, e.g., during rest and exercise conditions or 
exposure scenarios (4 12 . Such mechanistic descriptions of the observable exhalation 



kinetics give valuable insights into the relevance of the measured breath concentrations 
with respect to the endogenous situation and hence are mandatory to fully exploit the 
diagnostic potential of breath VOCs. 

Within this context, the main focus of this article will be on a modeling-based 
review of isothermal rebreathing, which has been proposed as an experimental tech- 



nique for estimating the alveolar levels of hydrophilic exhaled trace gases 13 15 . This 
class of compounds has been demonstrated to significantly interact with the water-like 
mucus membrane lining the conductive airways, an effect which has become known as 
wash-in/wash- out behavior. For further details we refer to [5 16 . As a phenomenologi- 
cal consequence, exhaled breath concentrations of such highly water soluble substances 
tend to be diminished on their way up from the deeper respiratory tract to the airway 
opening. The resulting discrepancies between the "true" alveolar and the measured 
breath concentration can be substantial (even if breath samples are drawn in a strictly 
standardized manner employing, e.g., CO2- or flow-controlled sampling) and will depend 
on a variety of factors, including airway temperature profiles and airway perfusion as 



well as breathing patterns 17 . 



In particular, the above-mentioned effect considerably departs from the classical 



Farhi description of pulmonary inert gas exchange 18 , on the basis of which the op- 
erational dogma has been established that end-tidal air will reflect the alveolar level 
and that arterial concentrations can be assessed by simply multiplying this value with 
the blood:gas partition coefficient Ab :a ir at body temperature. This "common knowl- 
edge" has first been put into question in the field of breath alcohol testing, revealing 
observable blood-breath concentration ratios of ethanol during tidal breathing that are 



unexpectedly high compared to the partition coefficient derived in vitro 13 . Similarly, 
excretion data (i.e., the ratios between steady state partial pressures in expired air and 
mixed venous blood) of highly water soluble compounds (including the MIGET test 
gas acetone) have been shown to underestimate the values anticipated by treating the 
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airways as an inert tube 19 20 . 



In a previously published mathematical model for the breath gas dynamics of highly 
soluble blood-borne trace gases, airway gas exchange is taken into account by separating 
the lungs into a bronchial and alveolar compartment, interacting via a diffusion barrier 
mimicking pre- and post-alveolar uptake M. This formulation has proven its ability 
to reliably capture the end-tidal breath profiles as well as the systemic dynamics of 
acetone in a variety of experimental situations and will be used here for illuminating 
the physiological processes underlying isothermal rebreathing tests as carried out in the 
literature. For comparative reasons, the illustration in the sequel will mainly be limited 
to acetone, with possible extrapolations to other highly soluble VOCs indicated where 
appropriate. 

2. Methods 

2.1. Experiments 

Extensive details regarding our experimental setup are given elsewhere (9,21,22 . Here, 
we will only briefly discuss the parts relevant in the context of isothermal rebreathing. 

The rebreathing system itself consists of a Tedlar bag (SKC, Dorset, UK) with 
a volume of Vb ag = 3 1, that can directly be connected to a spirometer face mask 
covering mouth and nose. From the latter, end-tidal exhalation segments are drawn 
into a Proton Transfer Reaction Mass Spectrometer (PTR-MS; Ionicon Analytik GmbH, 
Innsbruck, Austria), which allows for VOC detection and quantification on a breath- 



by-breath resolution as described in 21 . The rebreathing bag is warmed to 37 ± 1°C 



by means of a specially designed outer heating bag (Infroheat Ltd., Wolverhampton, 



UK), cf. 15 . Heating is intended to assist the thermal equilibration between the 
alveolar tract and the upper airways and prevents condensation and subsequent losses 
of hydrophilic VOCs depositing onto water droplets forming on the surface wall of the 
bag. End-tidal acetone concentrations are determined by monitoring the protonated 
compound at m/z = 59 (dwell time 200 ms). Additionally, we routinely measure the 
mass-to-charge ratios m/z = 21 (isotopologue of the primary hydronium ions used for 
normalization; dwell time 500 ms), m/z = 37 (first monohydrate cluster H30 + (H20) 
for estimating sample humidity; dwell time 2 ms), m/z = 69 (protonated isoprene; 
dwell time 200 ms), m/z = 33 (protonated methanol; dwell time 200 ms) as well as 
the parasitic precursor ions NH4 and at m/z = 18 and m/z = 32, respectively, 



with dwell times of 10 ms each. In particular, following 15 , the pseudo concentrations 
associated with m/z = 32 determined according to standard PTR-MS practice (see, 
e.g., [23| Eq. (1)]) will be used as a surrogate for the end-tidal oxygen partial pressure 
Po 2 (relative to an assumed nominal level at rest of 100 mmHg). Similarly, calibrated 
pseudo concentrations corresponding to m/z = 37 are considered as an indicator for 
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absolute sample humidity C wat er, see [9] for further details. Partial pressures Pco 2 °f 
carbon dioxide are obtained via a separate sensor. Table [T] summarizes the measured 
quantities used in this paper. In general, breath concentrations will always refer to 
end-tidal levels. 



Variable 


Symbol 


Nominal value (units) 


Acetone concentration 


^measured 


1 (Mg/1) 


24 




CO2 partial pressure 


Pco 2 


40 (mmHg) 


25 




Water content 


Cwater 


4.7 (%) 


2G 




O2 partial pressure 


Po 2 


100 (mmHg) 


25 





Table 1. Summary of measured breath parameters together with some nominal values 
in end-exhaled air during tidal breathing at rest. 



In total, five normal healthy male volunteers (age range 26-55 years, median 35 
years) were recruited and gave written consent to participate in a single cycle rebreathing 



protocol as described in section |3.2| All phenomenological results were obtained in 
conformity with the Declaration of Helsinki and with approval by the Ethics Commission 
of Innsbruck Medical University. 



2.2. Physiological model 

A schematic sketch of the model structure is presented in figure [T] For the associated 
compartmental mass balance equations and the underlying model assumptions we refer 



to section Appendix A as well as to the original publication [9]. Briefly, the body 
is divided into four distinct functional units, for which the underlying concentration 
dynamics of the VOC under scrutiny will be taken into account: bronchial/mucosal 
compartment (Cb ro ; gas exchange), alveolar /end-capillary compartment (Ca; gas 
exchange), liver (Cu v ; endogenous production and metabolism) and tissue (C^; storage). 
The nomenclature is detailed in table IAU 

The measurement process is described by 

Cmeasurcd Cbro; (1) 

i.e., we assume that the measured (end-tidal) VOC concentration reflects the bronchial 
level. Dashed boundaries indicate a diffusion equilibrium, described by the appropriate 



partition coefficients A, e.g., Ab :a i r . The nomenclature is detailed in table |A1 
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Figure 1. Sketch of the model structure used for capturing dynamic VOC 
concentrations (C). Subscripts connote as follows: b&g-rebreathing bag; I-inhaled; 
bro-bronchial; umc-mucosal; A-alveolar; c'- end- capillary; liv- liver; tis— tissue; b— 
blood. The conductance parameter D reflects the net gas exchange between the 
bronchial and the alveolar tract over a complete respiratory cycle as described in 
the text. Note that D is close to zero during free tidal breathing at rest (i.e., the 
measurable breath acetone levels are almost exclusively determined by the respective 
levels in the conducting airways) and varies in parallel with ventilatory flow, cf. (JHJ. 
In particular, higher values of D correspond to an increased contribution of the alveoli 
to total pulmonary gas exchange. 



In particular, during tidal breathing at rest the present model replaces the 
familiar Farhi equation describing the steady state relationship between inhaled 
(ambient) concentration Ct, measured breath concentrations C^ easured , mixed venous 
concentrations C v and arterial concentrations C a during tidal breathing (cf. |18|), 

&a + cv c 

^measured — U A — — , , \*) 

Ab:air + g A ^ir 

with the expression (cf. [9, Eq. (21)]) 

r _ r r hro Ci + (1 - gbro)Cy r hro Ci + C a 

^measured — ^bro — x . ~~ (rfi\\ , ' W 

(1 - gbroM-/jA b :air + ?"bro ^(T)A b :air + ?"bro 

Here, g bro ~ 0.01 is an estimate of the effective fractional bronchial perfusion (which 
might be substance-specific), while r bro = Va/ '(^broOc) denotes the associated bronchial 
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ventilation-perfusion ratio. Note that ([3]) is only valid for highly water- and blood- 
soluble compounds such as acetone. 

The term -2(T)Ab :a ir reflects an explicit temperature dependence of airway gas 
exchange that has been incorporated into the model. More specifically, the factor 
z captures the change of VOC solubility in the airway mucosa and bronchial blood 
in response to fluctuations of the characteristic mean airway and bronchial blood 
temperature T (in °C). In j9j Eqs. (3)-(5)] the latter has been shown to be estimable by 
virtue of the measured sample humidity C water and typically ranges around 34° C during 
tidal breathing of room air at rest. It is assumed that the temperature dependence of 
Ab:air is proportional to the temperature change of the mucosa: air partition coefficient 
A muC :air over the temperature range considered, i.e. z is given by 

z(T) '■— A muc:a i r (T)/A muc:a i r (37°C). (4) 

Note that z decreases with increasing temperature and takes values greater than one 
for T below body temperature. 

A particular convenience of the above representation stems from the fact that, as 
a first approximation, the mucosa layer can be assumed to inherit the physico-chemical 
properties of water. Correspondingly, A muc:a i r can be estimated via the respective wa- 
ter:air partition coefficient, which is usually available from the literature (see, e.g., the 
extensive compendium in (27j). 

A direct implication of the wash-in/wash-out mechanisms discussed in the 
introduction is that during free tidal breathing at rest the amount of hydrophilic VOCs 
reaching the alveolar region with each inhalation will be approximately equal to the 
amount extracted with each exhalation. In other words, under these circumstances the 
bronchio- alveolar net exchange of such compounds over one breath cycle will be close 
to zero and a major part of the overall pulmonary gas exchange (e.g., up to 96% as in 
the case of acetone (16)) can be ascribed to the conducting airways. As a consequence, 
a concentration gradient develops along the respiratory tree, with higher levels of the 
respective VOC found in the proximal (alveolar) parts. Increasing ventilatory flow 
rates and volume minimizes equilibration times with the more distal mucus layers and 
hence counteracts the saturation and depletion of the air stream during inhalation and 
exhalation, respectively. Consequently, the location of gas exchange is continuously 
shifted towards the alveolar region and the contribution of the latter becomes more 
pronounced. In particular, the end-exhaled (measured) breath concentration approaches 
values closer to the alveolar level. 

Formally, this rationale can be incorporated into the model by including a diffusion 
term of the form _D(Cb ro — Ca) into the underlying mass balance equations, cf. [9]. 
This expression reflects the above-mentioned net gas exchange between the alveolar and 
the bronchial tract over a complete breath cycle. From the previous paragraph the 
conductance parameter D will be close to zero during free tidal breathing at rest and 
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increases with ventilatory flow, viz., 

D := k diGA max{0, V T - V T est } + fc diff , 2 max{0, V A - V^}, (5) 
where Va and Vt denote alveolar ventilation and tidal volume, respectively. The positive 



coefficients /cdiff,i can be estimated from experimental data, see [9] and table Al 

A primary advantage of the above treatment is that much of the temporal and 
spatial heterogeneity associated with highly soluble gas exchange can be captured 
within a unidirectional rather than a bidirectional flow framework as in |4j[5]. From 
an operational view, this drastically reduces the number of required model equations, 
thereby yielding parsimonious parameterizations as well as tractable model structures. 

3. Results and discussion 

3.1. Heuristic considerations 

As has already been indicated in the introduction, ^ is inappropriate for capturing 
experimentally obtained arterial blood-breath concentration ratios (BBR) of highly 
water soluble trace gases during free breathing at rest (i.e., assuming Cj = 0). For 
instance, in the specific case of acetone, multiplying the proposed end-tidal population 
mean of approximately 1 fig/l [24j with a blood:gas partition coefficient of Ab ;a ir = 
340 1 16 1 at body temperature appears to grossly underestimate arterial blood levels 



spreading around 1 mg/1 [28,29 . In contrast, pi) asserts that the observable arterial 



blood-breath ratio in this case is 

BBR = ° & = z(T)A b:air + r bro , (6) 

^ measured 

and will thus depend on airway temperature and airway blood flow as mentioned in 
the previous section. In particular, note that the BBR will be greater than the in vitro 
blood:gas partition coefficient Ab :a ir for T below body temperature, cf. Q. 

From (|6j) it is clear that the more soluble a VOC under scrutiny, the more drastically 
its observable BBR will be affected by the current airway temperature, with an inverse 
relation between these two quantities. This deduction is consistent with measurements 
by [ll] conducted in the field of alcohol breath tests, reporting a monotonous decrease 
of ethanol BBR values with increasing exhaled breath temperature. In the same 
contribution, BBRs of ethanol during normal tidal breathing were shown to typically 
exceed a value of 2500, which differs from the expected value Ab :a ir = 1756 by more than 
40%. For comparison, based on a mean characteristic airway temperature of about 



T = 34°C 30 and the water:air partition coefficient A muc:air provided by 27 p. 568], by 



substituting the parameter values in table Al we find that (|6| predicts an experimentally 



observable blood-breath ratio of ethanol in the range of 2250. Equation ^ also suggests 
that in the case of VOCs with an extremely high affinity for both blood and water (e.g., 
ethanol, methanol), the discrepancies between in vivo and in vitro BBR values might 
largely be removed by warming the airways to body temperature before the breath 
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sampling procedure. Contrarily, blood-breath ratios of less soluble VOCs, e.g., acetone, 
will additionally be affected by the comparatively large value of r bro . 

Apart from providing some experimental evidence for the validity of (JgJ) , these ad 
hoc calculations suggest that the common practice of multiplying the measured breath 
concentration C measure d with Ab :a ir to obtain arterial concentrations for highly soluble 
trace gases will result in an estimation that might drastically differ from the true blood 
level. 

Isothermal rebreathing has been put forward as a valuable method for removing 
the above-mentioned discrepancies. The heuristic intention leading to isothermal 
rebreathing is to create an experimental situation where the breath levels of highly 
soluble VOCs are not altered during exhalation due to loss of such substances to the 
cooler mucus layer of the airways. This can be accomplished by "closing the respiratory 
loop", i.e., by continuous re- inspiration and -expiration of a fixed mass of aiifjjfrom a 
rebreathing receptacle (e.g., a Tedlar bag), causing the airstream to equilibrate with 



the mucosa linings over the entire respiratory cycle 14 16 . Additionally, warming the 
rebreathing volume to body temperature will ensure a similar solubility of these VOCs 
in both regions, alveoli and airways. 

Formally, a model capturing the experimental situation during isothermal rebreath- 
ing can simply be derived by augmenting the model equations in section |Appendix A 



with an additional compartment representing the rebreathing receptacle, i.e., 

'Hag = VA(Cbro — Cbag) (7) 



dCbag T ~, 



dt 



and by setting C\ = Cb ag in (A.l). In the sequel, we will furthermore assume 
that the (mixed) venous blood concentrations stay roughly constant throughout the 
rebreathing periods investigated in the present context. This can be seen as a necessary 
physiological requirement backed by the fact that both blood and most body tissues have 
an enormous storage capacity with respect to hydrophilic compounds, implying that the 
corresponding concentrations within these compartments will vary extremely slowly in 
response to a change in inhaled concentration levels. Consequently, we may conclude 
that if the conducting airways are warmed to body temperature during rebreathing 
(i.e., if z(T) in Q approaches unity as temperature increases) the compartmental 



concentrations in the pulmonary tract will tend to a new quasi steady statd^] obeying 

^-yrebr /^yrebr /-nchr /^yrebr a /o\ 

^measured ^bro ^bag \ \ ■ \°/ 

-^biair ^biair 

% For the sake of completeness, it should be noted that a small part of the air within the rebreathing 



circuit outlined in section 2.1 is actually consumed by the PTR-MS device (roughly 40 ml/min), which, 
however, can safely be neglected in the present setting. 

§ We stress the fact that this is not a steady state in the strict mathematical sense, i.e., corresponding 



to a situation where all time-derivatives of the mass balance equations in section Appendix A 



exactly zero. Effectively, in a hypothetical situation where rebreathing could be continued for several 
hours, breath acetone concentrations would still vary in parallel with a slow increase of body tissue and 
mixed venous blood levels. 
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This is a simple consequence of substituting C\ with Cb ro in ^ according to the quasi 
equilibrium relation associated with ([7]). Similarly, the classical Farhi model yields 



/^(F.rebr 



fiv ,rebr fiv ,rebr ^-yf ,rebr /q\ 

^measured ~~ bag ~~ — ~\ ~~ 7 ' ^ ' 



yF,rebr ^F^ebr ^F^ebr 

^b:air ^b:air 

It is important to realize that if such quasi steady state conditions hold they will de- 
pend solely on the blood:air partition coefficient Ab :a ir at 37°C, thus rendering isothermal 
rebreathing as an extremely stable technique for providing a reproducible coupling be- 
tween VOC levels in breath and blood. Particularly, it theoretically avoids the additional 
measurement of ventilation- and perfusion-related variables that would otherwise affect 
this relationship, thereby significantly simplifying the required technical setup for breath 
sampling. However, as will be illustrated in the following, the major practical obstacle 
is to guarantee that a quasi steady state as in (fsb is effectively attained. 



For the purpose of comparing the qualitative implications of (p]) and ^ it 
is instructive to note that the corresponding ratios between measured rebreathing 
concentrations and end-exhaled concentrations during free tidal breathing at rest will 
follow an entirely different trend. To this end, note that while from ^ and by setting 
Ci = in g we find that 

^measured _^ CyMbiair __ A b:a ir + Va/ Qc /.^ 

^measured Cy/ (^b:air + Va/ Qc) ^b:air 

an analogous computation using (|8j) and ^ shows that 

^measured ^ (j; ~ QbTo)z(T) Ab :a ir + 7"bro 
Cmeasurcd (1 (?bro) Ab:air 

where z(T) > 1, cf. Q. For highly soluble trace gases, this observation constitutes a 
simple test for assessing the adequacy of the Farhi formulation regarding its ability to 
describe the corresponding exhalation kinetics. Indeed, for sufficiently large Ab :a ir, the 



right-hand side of (10) will be close to one, while the right-hand side of (11) suggests 
that rebreathing will increase the associated end-tidal breath concentrations. In other 
words, for this class of compounds a markedly non-constant behavior during the initial 
isothermal rebreathing period indicates that ^ will fail to capture some fundamental 
characteristics of pulmonary excretion. Such tests are of particular importance in the 
context of endogenous MIGET methodology (Multiple Inert Gas Elimination Technique, 
based on endogenous rather than externally administered VOCs |20|), as they might be 
used for detecting deviations of the employed gases from the underlying Farhi description 
(see also [22]). 

3.2. Single cycle rebreathing 



In this section we will discuss some simulations and preliminary experiments conducted 
in order to study the predictive value of isothermal rebreathing within a realistic setting. 
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For this purpose we will first mimic isothermal rebreathing as it has been carried out 



by various investigators 13, 14 



More precisely, rebreathing was instituted after t s tart = 2 min of free tidal breathing 
at rest by inhaling to total lung capacity and exhaling until the bag was filled, thereby 
providing an initial bag concentration which can be assumed to resemble the normal 
end-exhaled level, i.e., 

C bag (0) = C bro (0). (12) 

Rebreathing was then continued until either the individual breathing limit was reached 
or the CO2 partial pressure increased above 55 mmHg. The volume of the rebreathing 
bag was Vb ag = 3 1 according to section 



2.1 



Typical profiles of end-tidal acetone, water, CO2 and oxygen content during normal 
breathing and isothermal rebreathing at rest are shown in figure [2] These representative 
data correspond to one single healthy male volunteer from the study cohort described 



in section 2.1 An overview of the experimental outcome for all volunteers is given 
in figure [3] 

As has been explained above, Po 2 can be derived by scaling the end-tidal level of 
the PTR-MS pseudo concentration signal at m/z = 32 to a basal value of 100 mmHg 
during free tidal breathing. Furthermore, breathing frequency / and tidal volume Vr 
may be simulated on the basis of the monitored values for Pco 2 an d Po 2 ■ Under iso-oxic 
conditions, after a certain threshold value is exceeded, both frequency and tidal volume 
are known to increase linearly with alveolar Pco 2 - Moreover, the corresponding slopes 
(reflecting the chemoreflex sensitivity of breathing) are dependent on the current alveolar 



oxygen partial pressure 25 , 31 . More specifically, hypoxia during rebreathing enhances 
chemoreflex sensitivity, yielding a hyperbolic relation between the mentioned slopes 
and Pq 2 . These findings result in a simple model capturing the chemoreflex control of 



breathing in humans 32 , which has been re-implemented in order to compute / and Vr 
from basal values during free breathing, cf. figure [2] Here, for simplicity it is assumed 
that end-tidal Pco 2 levels reflect those of the peripheral and central chemoreceptor 
environment. The ventilation rate follows from 

V A = f(V T -V D ), (13) 

using an instrument (mask) deadspace volume of Vd — 0.1 1. 

Since cardiac output can be expected to stay relatively constant during the 



rebreathing phase 15 , we fix its value at a nominal level of 6 1/min. Using the 



parameter values in table Al for a male of height 180 cm and weight 70 kg this completes 
the necessary data for simulating the above-mentioned rebreathing experiment. More 
specifically, the model response in the first panel of figure [2] is obtained by integrating 
the differential equations (A.1)-(A.4) and ([7|, setting C\ = Cb ag for t 6 [t s tart, ^end] 



and G\ = otherwise. The initial concentrations for each compartment as well as the 
underlying endogenous acetone production rate k pr = 0.055 mg/min for this specific 
experiment are derived from the algebraic steady state conditions at t — 0. 
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Figure 2. Representative outcome of an isothermal rebreathing experiment during 
rest. Data correspond to one single normal healthy male volunteer from the study 
cohort in section [2TT] (Volunteer 1, see also figure [3]). Isothermal rebreathing starts 
at istart = 2 min and ends at approximately t cn d = 4.6 min. Measured or derived 
quantities according to the experimental setup are shown in the first and third panel 
(red tracings), while data in the second and fourth panel correspond to simulated 
variables as described in the text (black tracings). 



From figure [2j the proposed model is found to faithfully reproduce the observed 
data. Further improvement of the goodness-of-fit might be achieved by employing 
parameter estimation techniques as described, e.g., in [8j[9] which, however, would be 
beyond the scope of this paper. 



In contrast, as can be anticipated from (10), when using the classical Farhi 
formulation for gas exchange in the the pulmonary tract (i.e., when treating the 
conducting airways as an inert tube) the given breath acetone profile cannot be 
capturec|jj] 

In particular, the presented data appear to consolidate the heuristic considerations 



in section 3^ and confirm that the alveolar concentration of acetone during free tidal 
breathing can differ from the associated bronchial (i.e., measured end-exhaled) level 



by a factor of more than 1.5. As has been explained in section 2.2, this is due to an 



|| Considering the fact that the Farhi formulation is included in the present model as a limiting case 
for (7bro = and D — > oo |9j Fig. 4], its associated output can be computed in a similar manner as 
described above. 
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effective concentration gradient between the conducting airways and the alveolar space. 
During isothermal rebreathing, this gradient slowly vanishes and causes the measured 
breath concentration to approach the underlying alveolar concentration. It should be 
noted that the latter as well as the simulated levels in liver and storage tissue stay 
almost constant throughout this period, which is consistent with the quasi steady state 
assumption made in (J8|. We also stress the fact that in order to simulate a similar breath 
profile by using the conventional Farhi formulation (i.e., by treating the airways as an 
inert tube), one essentially would have to postulate a temporarily increased endogenous 
acetone production during rebreathing, which, however, lacks physiological plausibility. 

For comparative reasons, in figure [2] and figure [3] we also display the simultaneously 
obtained PTR-MS concentration profiles of breath methanol, scaled to match the initial 
levels of breath acetone. Taking into account a methanol blood:gas partition coefficient 
of Ab:air = 2590 at body temperature [2], from (11) it can be deduced that for this 
compound the differences between concentrations extracted during free breathing and 
rebreathing primarily stem from the thermal equilibration between airways and alveolar 
tract. The corresponding rise in temperature is mirrored by a steady increase of sample 
water vapor C wa t e r, approaching an alveolar level of about 6.2%. In particular, the 
presented profiles for methanol show the necessity of including an explicit temperature 
dependence in models describing the rebreathing behavior of highly soluble VOCs. 
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Figure 3. Breath acetone/methanol profiles during isothermal rebreathing for all hve 
volunteers in the study cohort. Breath methanol profiles were scaled to match the 
initial levels of breath acetone. Rebreathing starts at t star t = 2 min. As soon as the 
individual breathing limit is reached (roughly between 4 and 5 min), the volunteers 
disconnect the rebreathing bag and continue with free tidal breathing at rest. 



As a last point in this section, we would like to emphasize the fact that all model 
parameters and initial conditions Pi governing the simulation in figure [2] can be ranked 
with respect to their impact on the observable breath acetone concentration during 
the isothermal rebreathing period [t s tart, ^end] by numerically approximating the squared 
L 2 -norm of the corresponding normalized sensitivities, viz., 

^cnd 2 

f \ f ( ^^measured(^) Pi \ , , 

5(P,) := J { dp, max.lyWlJ *' (M) 

^start 

Adopting the nomenclature in section Appendix A these indices were calculated for 
all effective compartment volumes, partition coefficients, initial concentrations, and 
for pi G {/c pr , ^met, 9bro, fcdiff,i, ^difi,2}- Among these, the decisive quantities are the 
blood:gas partition coefficient (negative influence; ?(Ab :a i r ) = 2.07) as well as the initial 
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concentrations and partition coefficients for the liver and tissue compartment (positive 
influence; s(C liy (0)) = 0.06; g(A b:liv ) = 0.06; <r(C tis (0)) = 1.39; ?(A b:tis ) = 1.39). In 
particular, the latter essentially define the mixed venous blood acetone concentration 



CV according to (A. 5). As this concentration stays almost constant throughout the 
entire experiment, the above sensitivity analysis results are in direct agreement with ([8]), 
predicting a rebreathing plateau that exclusively depends on Ab ;a ir and Cy. All other 
sensitivity indices take values below 0.015. 

3.3. Sequential rebreathing 

On the basis of the results from the previous subsection, in the following we will briefly 



discuss a sequential rebreathing protocol developed by 15 1. This regime aims at improv- 
ing the patient compliance of conventional rebreathing by repeatedly providing cycles of 
five rebreaths (postulated to last approximately 0.5 min) with intermediate periods of 
free tidal breathing lasting approximately 10 min. According to above-mentioned pro- 
tocol, isothermal rebreathing is again instituted by inhaling to total lung capacity and 
exhaling to residual volume into a Tedlar bag with a volume of T4>ag = 3 1. After each 
rebreathing cycle, the bag is closed, a small amount of bag air (< 100 ml) is measured 
and the volunteer starts the next cycle by exhaling to residual volume and inhaling from 
the bag. 

From the data in the last two panels of figure [2] it can be inferred that all 
physiological input variables will have returned to pre-rebreathing values within the 
10 min breaks. For simulation purposes, it will hence be assumed that their behavior 
during repeated rebreathing segments is identical to the profiles within the first 0.5 min 
of single cycle rebreathing. Values for the initial compartment concentrations as well as 
for the additional parameters are adopted from the previous subsection. These premises 
allow for simulating the evolution of breath acetone within a repeated rebreathing 
regime as displayed in figure |4| Here, the initial bag concentration at the onset of 
each individual rebreathing cycle is determined by the final bag concentration after the 
preceding rebreathing cycle, i.e., no fresh room air enters the bag. 

The bag concentration profile in figure [4] qualitatively resembles the data presented 



by 15 . However, what emerges from this modeling-based analysis is that in spite 
of steadily increasing bag concentrations (finally reaching a plateau level), the latter 
might not necessarily approach the underlying alveolar concentration as in the case 
of single cycle rebreathing. The major reason for this is a lack of complete thermal 
and diffusional equilibration between the airways and the alveolar region within the 
individual rebreathing segments. One potential way to circumvent this issue would 
be to reduce the desaturation and cooling of the airway tissues between consecutive 
rebreathing segments by keeping the intermediate time interval of free tidal breathing 
as short as possible (while simultaneously maintaining a regime allowing for comfortable 
breathing) . 
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The second panel in figure [4] displays the evolution of the predicted blood-bag 
concentration ratios during the course of experimentation. Note that the in vitro 
blood:gas partition coefficient Ab :a ir = 340 is never attained. This observation can offer 
some explanation for the discrepancies that continue to exist with regard to theoretical 



and experimentally measured ratios between blood and (rebreathed) breath levels 33 



Furthermore, the final plateau value and the observable BBR of acetone will vary with 
temperature (cf. (J6l)), which is consistent with similar observations made in the case of 



breath ethanol measurements 14 
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Figure 4. Simulation of a sequential rebreathing protocol according to 15 with 



intermediate pauses of 10 min characterized by free tidal breathing. The underlying 
model parameters correspond to the individual fit in figure [2] Dash-dotted lines 
represent upper and lower bounds for the bag concentration as well as for the observable 
blood-bag ratio with respect to changes in airway temperature T. These bounds were 
obtained by assuming that the airway temperature either instantaneously rises to body 
core temperature during the individual rebreathing periods or remains constant (i.e., 
at its initial level To) throughout the entire experiment. 



4. Conclusions 

Here we have successfully applied a previously published compartment model for the 
exhalation kinetics of highly soluble, blood-borne VOCs to the experimental framework 
of isothermal rebreathing. The proposed model has proven sufficiently flexible for cap- 
turing the corresponding end-tidal breath dynamics of acetone, which can be viewed 
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as a prototypical test compound within this context. Moreover, it allows for a quan- 
titative assessment of the discrepancies between in vitro and measured blood-breath 
ratios of such VOCs during the cyclical steady state observed in normal tidal breathing 
and can serve as a relatively simple mechanistic basis for describing some fundamental 
physiological processes associated with the observable VOC behavior in response to the 
above-mentioned experimental regime. 



Several practical implications emerge from this modeling-based analysis. Firstly, it 
is demonstrated that the classical Farhi setting will fail to reproduce the experimentally 
measured acetone exhalation data during isothermal rebreathing if a constant 
endogenous production and metabolism rate is postulated. This is due to the fact that 
airway gas exchange, being a major determinant affecting the breath levels of highly 
soluble VOCs, is not taken into account within this formalism. Furthermore, multiplying 
end-tidal breath concentrations during free tidal breathing with the substance-specific 
blood:gas partition coefficient Ab ;a ir will generally underestimate the true arterial blood 
concentration for highly blood and water soluble VOCs. 

Excessive hypoxia and hypercapnia are the main factors limiting the duration 
of the rebreathing maneuver and preventing a complete equilibration of VOC partial 
pressures in the alveoli and the conducting airways. From an operational point of view, 
our data indicate that even if isothermal rebreathing is continued until the individual 
breathing limit is reached, a quasi steady state according to ^ might not necessarily 
be attained. As can be deduced from figure [2]in the case of acetone, end-exhaled breath 
(or bag) concentrations extracted after about 0.5 min of rebreathing (corresponding to 
the common protocol of providing around five consecutive rebreaths) are still likely to 
underestimate the underlying alveolar level Ca- Analogous conclusions can be drawn 
for sequential rebreathing protocols designed to allow for a recovery of the volunteer 
during the individual rebreathing segments. 

A reliable extraction of meaningful breath levels for highly soluble VOCs by virtue of 
isothermal rebreathing hence appears to require more sophisticated setups incorporating 
the continuous removal of CO2 and replacement of metabolically consumed oxygen. 
Provided that the influence of chemical C02-absorption on the measured breath and 
bag concentrations is negligible, such setups might for instance be adapted from closed 
chamber techniques 34 or general closed-circuit anesthesia systems. 

Alternatively, alveolar concentrations might be extrapolated to some extent 
from partially equilibrated rebreathing samples, using numerical parameter estimation 
schemes for reconstructing VOC exhalation kinetics according to a given model 
structure. While further validation and data gathering needs to be carried out before 
such estimates can become practically relevant, the mechanistic descriptions discussed 
in this paper are intended as a first step towards achieving this goal. 
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Appendix A. Model equations and nomenclature 

This appendix serves to give a roughly self-contained outline of the model structure 
sketched in figure [l]|9|. The time evolution of VOC concentrations is captured by taking 
into account standard conservation of mass laws for the individual compartments. 

Local diffusion equilibria are assumed to hold at the air-tissue, tissue-blood and 
air-blood interfaces, the ratio of the corresponding concentrations being described 
by the appropriate partition coefficients A, e.g., Ab ;a ir- Unlike for low blood soluble 
compounds, the amount of highly soluble gas dissolved in the local blood volume 
of perfused compartments cannot generally be neglected, as it might significantly 
increase the corresponding compartmental capacities. This is particularly true for 
the airspace compartments. We hence use the effective compartment volumes Vb ro := 
Vbro + Vm UC A muc:a i r , Va := Va + K'A^air, V\ iv : = Vu v + Vii Vj bAb : ii v as well as V t i B := V t i s and 
neglect blood volumes only for the mucosal and tissue compartment. 

Values for the individual compartment volumes, temperature-dependent partition 
coefficients, and physiological variables such as cardiac output (Q c ) and alveolar 
ventilation (Va) are given in table Al , Fractional blood flows to the bronchial tract (quo) 



and the liver (qn v ) as well as the endogenous production (k pi ) and metabolic elimination 
rates (fc me t) are treated as constants in the context of isothermal rebreathing. 

According to figure [TJ the mass balance equation for the bronchial compartment 
reads 

dC*b 

^ Vbro = Va(Ci — Cbro)+-D(C , A — ^bro)+9broQc(Ca. — ^(T) AbiairCbro) , (A.l) 

where C\ denotes the inhaled (ambient/bag) gas concentration. Similarly, for the 
alveolar, liver and tissue compartment we find that 
dC A ~. 



and 



and 



dt 

dew .-, 
dt 

dC tis 
dt 



Va — D(C hro — Ca) + (1 — (Zbro)Qc(Cv — A^air^A.)) (A. 2) 

Vliv = kpr — &mctAb:livCliv + gn v (l — <7bro)Qc(Ca — Ab:li v Cliv), (A. 3) 

V tis = (1 - ftiv)(l - gbro)Qc(C a ~ A b:tis C tis ) , (A.4) 
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respectively. Here, 

Cv '■= <?livAb:livCliv + (1 — <?liv)Ab:tisCtis (A. 5) 

and 

Ca := (1 — <?bro)Ab:airCA + gbro^(^) Ab:airCbro (A. 6) 

are the associated concentrations in mixed venous and arterial blood, respectively. 
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Parameter 



Symbol Nominal value (units) 



Compartment volumes 

Bronchioles 

Mucosa 

Alveoli 

End-capillary 

Liver 

Blood liver 
Tissue 

Rebreathing bag 
Respiratory parameters 

Breathing frequency 

Tidal volume 

Alveolar ventilation 
Blood flows 

Cardiac output 

Fractional flow bronchioles 

Fractional flow liver 
Partition coefficients (37° C) 

Blood:air 

Mucosa: air 

Bloodtlivcr 

Blood:tissue 

Metabolic and diffusion constants 



A 



Vbro 
^muc 

V A 
V* 

V!iv,b 

V tis 

/ 

V T 
Va 

& 

9bro 
QlW 

Ab:air 
muc:air 
Ab:liv 
Ab:tis 



O.f (1)° 
0.005 (1)° 

4.1 (l) a 
0.15 (l) b 
0.0285 LBV (l) a 
1.1 (l) c 
0.7036 LBV (l) a 

3 (ir 

12.5 (tides/min)< i 
0.593 (\) d 
6.2 (l/min) e 

6 (1/min)^ 
0.01 s 
0.32 Q 

340^ 
392^' 
1.73 j 
1.38 h 



Linear metabolic rate 


^met 


0.0074 (l/kg°- 75 /min) 


Endogenous 


production 




0.19 (mg/min) 


Stratified conductance 


D 


(1/min) 


Constant ( 


5 






14.9 (min- 1 ) 


Constant ( 


5 




&diff,2 


0.76 



Table Al. Basic model parameters and nominal values during rest. LBV denotes the 
lean body volume in liters calculated according to LBV = —16.24 + 0.22 bh + 0.42 bw, 
with body height (bh) and weight (bw) given in cm and kg, respectively (6j; ° (6j; 



35 ; c 36 ; d 32 ; e cf. (1131) ; / 37 



25 



16 ; 1 27 
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